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Localized growth modes, dynamic textures, and upper critical dimension for the 
Kardar-Parisi-Zhang equation in the weak noise limit 
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A nonperturbative weak noise scheme is applied to the Kardar-Parisi-Zhang equation for a growing 
interface in all dimensions. It is shown that the growth morphology can be interpreted in terms 
of a dynamically evolving texture of localized growth modes with superimposed diffusive modes. 
Applying Derrick's theorem it is conjectured that the upper critical dimension is four. 
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There is a current interest in the general morphology 
and in particular the scaling properties of nonequilibrium 
models in statistical physics. There is, moreover, a need 
to develop general methods beyond perturbation theory, 
renormalization group theory, mode coupling theory, and 
numerical simulations, which permit an analysis of these 
often intractable problems. 

The purpose of this Letter is two-fold. On the one 
hand, I should like to draw the attention of the com- 
munity to the availability of a weak noise approach to 
stochastic equations driven by Gaussian noise which is 
based on a principle of least action and which allows 
a discussion of stochastic processes in terms of classical 
equations of motion and, moreover, provides the Arrhe- 
nius factor or weight of a specific kinetic transition; the 
method is summarized below. Secondly, I apply the weak 
noise approach to the Kardar-Parisi-Zhang (KPZ) equa- 
tion for the kinetic growth of an interface in arbitrary 
dimensions and show that the growth morphology can 
be interpreted as a dynamical network of growth modes 
with superimposed diffusive modes. As a corollary, us- 
ing Derrick's theorem, it is finally shown that the texture 
of growth modes does not persist above four dimensions, 
indicating that the upper critical dimension of the KPZ 
equation is d = 4. Details and further developments will 
be discussed elsewhere. 

The variationally-based weak noise method dates back 
to Onsager and has since reappeared as the Freidlin- 
Wentzel theory of large deviations H, |E ^ and as the 
weak noise saddle point approximation to the functional 
Martin-Siggia-Rose scheme 0, . The point of departure 
is the Langevin equations for the set of stochastic vari- 
ables Wn, 



1, ■ ■ ■ N, driven by white Gaussian noise 



dwp 

{rin'nm){t) = ASn„iS{t), 



2 ^ 



2 ^mn^jn^pn 



(1) 

(2) 



where Fn(wp) is the drift, Gnmiwp) is accounting for 
multiplicative noise, V„ = d/dwn and A is the explicit 



noise strength; sums are performed over repeated indices. 
The associated Fokker-Planck equation for the probabil- 
ity distribution P(w„,t) then has the form 0,f 



(3) 



where the symmetrical noise matrix Kpm{wp) 

Gpn [Wp ) Gmn [Wp ) . 

Introducing the WKB ansatz 

P{wn,t) cx exp[-S'(w„, i)/A], 



(4) 



the Hamilton- Jacobi equation dS/dt + H = 0, pn — WnS 
follows to leading order in A with Hamiltonian 



The Hamilton equations of motion are 
'^'"^ --F +K V 



dt 
dt 



^Prn^nF,n ^PmPq^^Rmq, 



(5) 

(6) 
(7) 



determining classical orbits on the energy surfaces given 
by _ff in a classical {wn,Pn) phase space. Finally, the 
action 5* is given by 



dtd^w Pn 



dw„ 
~dt 



- HT. 



(8) 



The weak noise scheme bears the same relationship to 
stochastic fluctuations as the WKB approximation in 
quantum mechanics, associating the phase of a wave 
function with the action of a classical orbit. In addi- 
tion to providing a classical orbit picture of stochastic 
fluctuations and thus allowing the use of dynamical sys- 
tem theory, the method also yields the Arrhenius fac- 
tor P cx exp(— 5'/A) for a kinetic transition to Wn in 
time T. Here the action S serves as the weight in the 
same manner as the energy E in the Boltzmann factor 



2 



P oc exp{—E/kT) for equilibrium processes. This com- 
pletes the brief review of the nonperturbative weak noise 
scheme; for details see e.g. Ref. |3| 

The KPZ equation is a field theoretic Langevin equa- 
tion descri bing the nonequilibrium growth of an interface 



— ^vV^h+ -Vh-Vh 
at 2 

(w)(r,t) = A<5''(r)^(t). 



F 



(9) 
(10) 



Here h{r,t) is the height of the growth profile, v a dif- 
fusion coefficient, A a growth coefficient, F an imposed 
drift, and rj a locally correlated white Gaussian noise 
of stren gth A. Dynamic renormalization group (DRG) 
studies |lll 113 . Iis ll indicates that the KPZ equation 
conforms to the dynamical scaling hypothesis with long 
time -long distance correlations {hh){r,t) = r^f$(t/r^), 
characterized by roughness exponent ^, dynamic expo- 
nent z and scaling function The KPZ equation is, 
moreover, invariant under the Galilean transformation, 
r ^ r - Xu°t, h^h + u°-r, F^F + (A/2)u° • u^, 
and the slope V/i = u ^ u + u" which, implying the 
scaling law C + ^ — 2, delimits the KPZ universal- 
ity class. In cZ = 1 the stationary distribution of h is 
P{h) cx exp[-(z//A) / dx (V/i)^] Q, yielding ( ^ 1/2 
and z = 2 — r/ = 3/2. In d > 2 the DRG implies a kinetic 
phase or roughness transition at a finite A from a smooth 
phase with z — 2, and C = (2 — d) /2, the linear Edwards- 
Wilkinson (EW) case for A = [la, to a rough phase 
with still debatable scaling exponents z and C = 2 — z, 
see e.g. Refs. ^8, 19, 20,^]. It has, moreover, been 
conjectured that d = 4 is an upper critical dimension 
beyond which the KPZ equation exhibits EW behavior 
|2^^, 24]. 

Applying the nonlinear Cole-Hopf transformation 
'w{r,t) — exp(A/i(r, i)/2zy) the KPZ equation (jS)) takes 
the form 



dw 2 ^ 

— = J/ V w- —wF + —wr], 
at 2v 2v 



(11) 



with multiplicative noise. In a moving frame w is changed 
according to w ^ w exp[(A/2i^)u° • r]. Note that in the 
noiseless case for 77 = ((TT|l reduces to the linear dif- 
fusion equation permitting a rather complete analysis of 
the deterministic KPZ equation for a relaxing interface 
[lllE^j l. The Cole-Hopf equation (|ll|l with multiplicative 
noise forms the basis for the mapping of the KPZ equa- 
tion to a model of directed polymers (DP) in a quenched 
random medium which by means of the replica method 
relates the roughness transition to a pinning transition 
in the DP model [MII^. 

In recent work I have applied the weak noise method 
to the KPZ equation or, equivalently, the noisy Burg- 
ers equation in d = 1 [9|, I23. Here the method read- 
ily yields a consistent Galilean invariant dynamical pic- 



ture of a growing interface in terms of propagating do- 
main walls with superimposed diffusive modes. The lo- 
calized domain walls account for the growth; the diffu- 
sive modes form a subdominant background governed by 
EW dynamics and not contributing to the growth. The 
method, moreover, identifies scaling exponents, univer- 
sality classes, and associates the dynamic exponent z 
with the domain wall dispersion law or the Hurst ex- 
ponent H = 1/z with the anomalous superdiffusion of 
growth modes, see also Ref. [l^ . 

Drawing on the insight gained in d = 1, the weak 
noise method is here _applied to the KPZ equation in ar- 
bitrary dimensions |2]|. Unlike the d = 1 case, where 
the local slope field u = V/i is the natural variable, the 
Cole-Hopf diffusive field w is more convenient for analy- 
sis in d > 1. Referring to the general weak noise scheme 
and making the assignment, w„ w, Pn p, Knm — *■ 
{X/2ufw^S'^{r - r'), and F„ ^ -2[iyV^w - {X/2iy)wF], 
the Hamiltonian density is given by 

7^ =]3[i/V„V„-i//c2]w + ifc2^V, (12) 

and the canonical field equations take the form 
dw 



— — — v\7^w — vk^w -\- k^w^p, 
at 



(13) 



^ = -i/V^p + vk^p - klp^w. (14) 

Here k — {\F /2v'^Y^'^ and fco — \/2v define two char- 
acteristic inverse length scales. Using the equation of 
motion H13|l the action is 



S{w,T) = -kl 



w,T 



d'^xdt 



2 2 
w p , 



(15) 



Wi,0 



and P{w,T) cx exp(— 5(u',T)/A) yields the transition 
probability. 

The action and as a result the equations of motion 
are invariant under the combined Galilei transformation 
w — + w exp[(A/2z/)u° ■ r], p ^ pexp[— (A/2i^)u'' • r]. The 
equations of motion determine orbits in the canonical 
{w,p) phase space from an initial configuration Wi to a 
final configuration w traversed in time T with the noise 
field p as a slaved variable. The orbits lie on the con- 
stant energy surfaces H = J d'^'x H. The action evalu- 
ated along the orbit then yields the transition probability 
P{wi — *■ w,T). The formalism is symplectic (canonical) 
and it is an easy task to perform canonical transforma- 
tions to the height field h or the local slope field u — V/i. 

The growth of the interface is due to the propagation of 
localized modes across the system. As in the d — 1 case 
the first task is thus to identify the relevant excitations 
and connect them in a dynamical network in order to ob- 
tain a consistent growth morphology. In the static limit 
the field equations H13|) and H14|) assume the symmetrical 
form: vV^w = vk'^w — k'^w'^p and i/V^p = vk'^p — k'^p^w. 
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On the noiseless manifold, p = 0, and the noisy mani- 
fold, p = vw, the static equations reduce to the linear 
diffusion equation and the nonlinear Schrodinger equa- 
tion (NLSE), respectively. 



V^w - '-2 



for p — Q , 
^ for p — vw. 



(16) 
(17) 



The diffusion equation (|15|) admits the radially sym- 
metric solution w+(r) oc r^~'^^'^Ii-d/2{kr), where I{z) is 
the modified Bessel function j^^. For small r, w+{r) 
est. For large r, w+{r) cx r'^^"'^)/^ exp(A:r), yield- 
ing the asymptotic height field h+{r) = (2j//A)((l — 
d)/2)\nr + kr) « {2v/\)kr^ and the outward-pointing 
vector slope field u+(r) = {2v / X)kY / r of constant mag- 
nitude (2v/\)k. In rf = 1, w+{x) oc cosh(fca:), giving 
rise to the height field h+{x) = (2t^/A) lncosh(A:a;), and 
the slope field u+{x) = Wh = (2i//A)A:tanh(A:a:), i.e., the 
right hand domain wall solution of the static Burgers 
equation i/V^u = —XuVu [s^, Since the solution 

w-i- lives on the noiseless manifold p = and the action 
S-^ = it carries no dynamics. 

The NLSE ()17(l admits a radially symmetric bound 
state w (r) falling off as w_(r) cx r^^"'')/^ exp(— fcr) for 
large r [32, 113, 0|- origin w_(0) = Udk/kg, 

where ad depends on the dimension. Numerically, oi = 
1.41(72), 02 = 2.21, 03 = 4.34. In the limit d ^ 4 
the amplitude diverges, the width vanishes, and the 
bound state disappears for d > 4. The asymptotic 
height field h^{r) = —{2v/X)kr and the inward-pointing 
slope field u_ = — (2i//A)fcr/r of constant magnitude 
{2v/X)k. In c? = 1 the NLSE admits the soliton solu- 
tion w^{x) = ai{k / k^) coshr^ {kr) yielding the height 
field h_{r) ~ — (2i^/A) Incosh(fcr) and slope field u_ = 
V/i_ — -~(2t^/A)fc tanh(fcr), i.e., the noise-induced left 
hand domain wall solution |35|, [Sg ■ Here the solution w- 
is associated with the noisy manifold p = uw and carries 
according to ifT^ the action S- = {kl/2)T J cfix wt- In 
Fig- n the radial solutions of the NLSE are depicted with 
parameter choice k = 1 and ko — 1 in d = 1,2, 3, 3.5. 

At large distances the slope fields u± associated with 
w± approach vector fields of constant magnitude 2{y/X)k 
and the boundary condition u = 0, corresponding to a 
flat interface, is implemented by combining a set of modes 
with appropriately chosen amplitudes /c. In a charge lan- 
guage, connecting modes with positive {k > 0) and neg- 
ative charges (fc < 0) and enforcing charge neutrality 
ki = 0, the boundary condition is automatically en- 
forced. 

The construction of the network is implemented in 
terms of the slope field u. Assigning a dilute network of 
static modes at positions , z = 1 • • • with charges ki the 
total slope field is u(r) = (2i//A) X), Vu;,(|r-rO|)/-u;i(|r- 
r"|), Xj ki — 0. In the vicinity of the mode position 
the slope field is shifted by u° = {2v/X) ^^^i Vwi(|r° — 
r"|)/-u;i(|r° — r?|) and the Z-th mode is assigned the ve- 




FIG. 1: The radially symmetric bound states of the NLSE 
are depicted for A: = feo = 1 in d = 1, d = 2, d = 3, d = 3.5. 



locity V; = — Au°. Using the asymptotic expressions for 
the modes and introducing a core radius e or order k~^ 
the self- consistent dilute dynamical network constituting 
the growth of an interface together with the associated 
action is given by 



h{v,t)^-Y,h^/{v-vm^+e, (18) 

i 

n{t) = f Mt')dt' + rl V k, = 0, (19) 
do 



V, (t) = -2u^ki 



r,{t)-ri{t) 
\Mt)~ri{t)y 



(20) 



^x wl (21) 



i,ki<0 



In Fig. 121 we depict a 3D snapshot of a 4-node growth 
morphology in d = 2. 

In the linear case for A:o = the equations of motion 
(|13|l and (|14|l admit extended diffusive modes with dis- 
persion uj = vk^ , corresponding to the EW universality 
class, as discussed in d = 1 case in Refs. |^ H^. The 
linear modes do not contribute to the growth and are su- 
perimposed on the dynamical network in the KPZ case. 

On the basis of mode coupling and DP theory [j^ 
1^ it has been argued that the upper critical dimension 
for the KPZ equation is dc — 4. In the present context 
dc is associated with the absence of growth modes above 
d = 4 as indicated by the numerics of the NLSE. Here I 
present a more rigorous argument for the existence of an 
upper critical dimension using Derrick's theorem 
based on constrained minimization. 

The NLSE (|17|l can be derived from the variation 
of the free energy F = K + {l/2)k^N — k^I, where 
K — (1/2) J d'^x (Vw)^ is the deformation energy, N = 



4 
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FIG. 2: 3D plot of a 4-node height profile with nodes at 
(±20, ±20) and charges 2.0, -1.5, -1.0, and 0.5 (in units 
of 2i^/\) 



J dfix the norm, and I = J d'^x w'^ the interaction, 
i.e., SF/Sw = yields V^w = k^w — k^w^. More- 



over, multiplying the NLSE by w and integrating over 
space yields the first identity: —2K = k'^N — kgl. Un- 
der the scale transformation w{r) w{iJ,r) one infers 
K fi'^-^K, N -> ^I'^N, and / fi'^I and subject 
to constrained minimization di^/d/i = the second 
identity: id-2)K + {k^ /2)dN - {kl/4)dl = 0. Eliminat- 
ing K from the identities one infers k'^N = {k^ / 4) {4— d) I . 
Since > and / > it follows that d < 4 in order for 
a bound state to exist with finite norm. 

In the present Letter I have summarized a general non- 
perturbative variationally-based weak noise approach to 
white noise driven stochastic processes. The method has, 
moreover, been applied to the KPZ equation in arbitrary 
dimensions yielding a Galilean invariant self-consistent 
dynamical network of modes accounting for the kinetic 
growth of the interface. In d = 1 the results agree with 
earlier findings, i.e, a network of matched domain walls 
whose dispersion yields the dynamic exponent z = 3/2; 
in d > 1, the detailed scaling properties remain to be 
worked out. Finally, based on the dynamical network 
representation and constrained minimization I have given 
an argument for d = 4 as the upper critical dimension for 
the KPZ equation. 
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